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^^ ' Abstract 

, -» ' We present a exactly soluble model for financial time series that mimics 

^— «, , the long range volatility correlations known to be present in financial data, 

r^ ' Although our model is 'monofractal' by construction, it shows apparent 

O^ I multiscaling as a result of a slow crossover phenomenon on finite time 

scales. Our results suggest that it might be hard to distinguish apparent 

C^ ' and true multifractal behavior in financial data. Our model also leads to 

S , a new family of stable laws for sums of correlated random variables. 

1-^ . Many time series exhibit interesting scaling properties. This means that 

r^ ' if x(t) denotes the time series, the probability distribution of the variations 

O ■ StX = x{t + T) — x{t), rescaled by a lag-dependent factor ^(T), can be written 

as: 



Vh ' where J-'{u) is a time independent scaling function. For example, if x{t) is 



constructed by summing independent identically distributed random variables 
with finite variance, one has (,{T) = o\/T and T[u) = exp(— u^/2)/\/27r for 
large T. Note that often the 'time' is actually a space coordinate, as it is the 
case in the analysis of turbulent velocity fields (where x is the fluid velocity) 
[y , or fracture surfaces (where x is the height of the profile) [g| . Equation (^ 
implies that all moments of 5tx that are finite, scale similarly: 

m, ^ {\STxn = A^arr, (2) 

where Ag is a q-dependent number. Very often, ^(T) behaves as a simple power- 
law: ^(T) ex T^. In this case of a monofractal process, one therefore has, 
m„ ex T^' , with („ = qd^. 



This is however not the only possibihty, and in some cases, one can observe 
multifractal scaling, in the sense that mq ex T^', with (q ^ q(. Such a possibihty 
has been advocated for turbulent velocity fields 0, 0, 0, 0] and, more recently, 
for financial time series |7[ ^, ||, |l^ . In the case of turbulence, there is strong 
theoretical evidence in favor of such a multifractal behavior pi, 0]. One can 
actually analytically derive a non trivial function ^q within a simple ('passive 
scalar') model, which is thought to retain some essential features of real tur- 
bulence ||ll[. The situation is much less clear in the case of financial markets, 
where the only evidence is based on the empirical analysis of the moments of 
several time series (typically currencies or stock indices) . The idea of multiplica- 
tive cascades, which is at the heart of the arguments in favor of multiscaling in 



turbulence, is not easily applicable to price time series (see, however, 12 ). 

In this note, we study to which extent empirical studies on multiscaling 
behavior in finance are sensitive to crossover behavior that results in apparent 
multiscaling, even though the studied process is a monofractal. To that end we 
present a soluble model that is based on the study of financial time series. In 
the model, the 'volatility' (or the variance) of the elementary price increments 
is a random variable with long range correlations, which have been shown to 
be present in financial time series |lj, |l^, |lj, |5| |l^. The model is an exact 
monofractal, but nevertheless it leads to an apparent multiscaling behavior p7| , 
Q. As we argue below, one finds effective exponents (q 7^ C? due to a very 
long crossover effect, which leads to a systematic negative correction to the 
true asymptotic behavior C^ = q(^. The correction grows with q and results 
thereby in a nontrivial functional form of (^q. The numerical simulation of such 
a model, which mimics quite well the observed behavior of real prices, accurately 
reproduces the published data in favor of multiscaling in financial markets. 

The model that we propose is also interesting in its own right. As a function 
of the strength of the correlations, we find a transition between a simple Gaus- 
sian behavior for the scaling function JT (together with the usual vT scaling for 
S,{T)) for weak correlations, to a new family of stable laws (with a non trivial 
scaling of ^(T)) for strong enough correlations. This adds to the very few cases 
where the limit distribution for sums of correlated random variables is exactly 
known. 

Our model is the following: we consider that the time series is built by 
summing random variables: 

^ f 

a;(t) = ^fefc N=-, (3) 

fc=i 

where t is a microscopic time scale. The elementary increments Sxk are assumed 
to be given by the product of two independent random variables, a 'sign' e^ 
and an amplitude ak'- Sxk = \crk\^k- The e^ are furthermore assumed to be 
independent Gaussian random variables of variance unity, af, therefore is the 
(random) variance of the elementary increments. We choose the ct/j's to be 



Gaussian random variables of zero mean p9| | , with a correlation function given 
by {(Jk(Jk+t) = C'd^D- AH moments niq of 5xk are therefore finite, and the even 
ones given by m^ = [{q - iy:]'^C''/'^{Q). 

The correlation function C(|£|) will be chosen to be a power law for large 
arguments: C(|£|) ~ ^\i\~^ . From several studies of financial markets, one 
knows that the variance of the price increments is indeed also random, with a 
very slowly decaying time correlation function: 



cr,- cr,- 



^)-{aty = 2C\\l\)^,^^-^, (4) 



2^ 



where the exponent v is found to be on the order ol v = 0.1 - 0.3 for different 
markets O, nj, HS, 16, 13]. The important point here is that v < 1/2. 



More precisely, we will use the following explicit representation of the cr's: 

N/2 . , ]^ 

-^ - ^ E wPi ^ (-™^"^ + <„e-"^) • (5) 



=1 

where the Zm's are independent complex Gaussian variables of unit variance. 
In the large iV limit, the resulting correlation function Cii) is well defined and 
decays, for large £, as a power-law with 7 = s^(y) cos{t:v/2) /it while C(0) 
tends to sgTT'^-Vz/ [|o|. 

We now turn to the calculation of the cumulants Cq of Pn{x) of x{t), as 
given by Eq. (0). We will show that these cumulants scale anomalously with N 
as soon as n > 2 (for < ^ < 1/2). 

After making a gauge transform e^ — > sign(cri)ei, one finds plj |: 

g-G«(^) ^ f dxPN{x)e-"'' (6) 

i=i ]=i j,k=i 



nfdejd(7-j\ 1 



The Gaussian integrals can be easily performed, and leads to the following 
expression for the characteristic function G{z): 

N/2 

Gn{z) = -Trlog(l + z^C) = Y. l°g(l + ^^C{m)), (7) 

rn— 1 

where the bold characters is used for matrices, and where C{m) are the eigen- 
values of the matrix C. From the very construction of the cr's, one finds that 
C{m) = SQ{2Trm/NY'^, each of which is twofold degenerate. 

Expanding Gn{z) in powers of z leads to the cumulants Cq{N) of Pn{x). All 
odd order cumulants are zero, while even order cumulants are given by: 

N/2 

Cq{N) = 2{q- 1)1 Y^Cimy/'. (8) 

m— 1 



Let us first analyze the case < i^ < 1/2. In the large N limit, the sum over m 
is convergent when g > 4 and leads to cumulants which do not scale as N: 

c,c.(£) sliq-l)lY^mi('^-^\ (9) 

^ ^ m=l 

while ioT q — 2 one finds exactly C2 = TrC = NC{Q). The normalized cumulants 
Cq/c2 therefore behave as N^'^"^^ for g > 4 and vanish for N -^ oo. This 
means that the distribution of x/^fN indeed tends to a Gaussian for large N . 
However, the approach to the Gaussian is slower than for sums of independent 
random variables. In particular, the kurtosis C4/C2 of Pn [x) decays anomalously, 
as N^"^^ , for Q < v < 1/2. For larger values of v (i.e. when the volatility 
correlations are weaker), one recovers the usual scaling of the kurtosis as \/N 
that holds for independent increments. Such an anomalous decay of the kurtosis 
with time was first reported for financial time series in pq , na| . 

The important outcome of the above calculation is that the moments vriq of 
the distribution Pn{x) are not simple power-laws, but sums of power laws with 
similar exponents. For example: 

7714 = AiflN^'^'' + A2^2N^, (10) 

7776 = ^6,0^'-^" + A4,2iV^"'" + A2,2,2^^^ (H) 

where the A's are some coefficients. If v is small, these sums of power-laws can be 
fitted on many decades with an effective exponent Cg such that ruq ex iV'»'' . The 
exponent C,q is less than q/2^ and more and more so as q increases. However, the 
true asymptotic behavior predicted by our model is C,q = 9/2- In fact, Eqs. ( p^ 
and (jlj) show that our model (and possibly also real financial data) is better 
characterized by the cumulants than by the moments. 

In order to illustrate this point numerically, we have generated a surrogate 
time series in a way closely related to the above model. Instead of writing 5xi = 
|(Ti|ei, we have chosen to take 5xi = exp(cri)ei. This leads to a more realistic 
time series as compared with real data from financial markets, without changing 
the crucial feature of the above model, i.e. the very slow decay of the volatility 
correlations. In particular, the distribution of the volatility has a positively 
skewed, log-normal shape. The length of our surrogate time series was taken to 
be comparable to those analyzed previously. The moments 771^ = (|(5Ta;|^) are 
plotted as a function of T for different q's, for the choice v = 0.2 (see Fig. 1). The 
interval of T was chosen to be r = 1 < T < 6000, again comparable to the region 
investigated in previous studies [0, ||, ||, 0. The power law fits are extremely 
good, and lead to a function Qq bending downwards as q increases, shown in 
Fig. 2. For our choice of parameters, the numerical values of Cg actually match 
precisely those reported in [1^. We have also checked that the same model, but 
without volatility correlations, leads very precisely to Cg — Cg = 9/2- 
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Figure 1 : Simulation of our surrogate financial time series and determination of 
the scaling of the moments as a function of window size T. From top to bottom 
(|(5Ta;|'^) for q=l - 10. Each moment scales relatively well with T as indicated 
by the 10 linear (log-log) fit lines. In this simulation the log- volatility follows a 
correlation Gaussian process as defined in the text with ly = 0.2. The variance of 
this process has been set so that the kurtosis at the unit time scale be kq — 65. 
This is a typical value for high frequency financial data. The simulated data set 
contained 500 000 points. 




Figure 2: Slope C,q of the fitted data of Fig. 1. The deviations from the true 



asymptotic scahng (q ~ q/2 are quite clear, 
those reported in px[ |. 



The value of C^ precisely match 



We now turn to the case i^ < 0. This corresponds to a non stationary process 
for the volatility, which typically grows with i as i~^l'^ . More precisely, from 
Eq. (|^), one can show that ((crfe — CTfc+f)^) oc l"^ . In this case, after changing 
variables to x = xyK^^ and z — z\J(N JI-kY^'^ ^ one finds that the asymptotic 
distribution of x has a characteristic function given by: 



G,(z)=^log(l 



~2 h 

z m 



(12) 



(we have set Sq = Ij which amounts to a change of scale in x.). 



The above result means that after rescaling by a factor VN^~'^, the sum of 
(strongly) correlated random variables converges to a non-Gaussian distribution 
Tu{x), obtained as the Fourier transform of the exponential of —G^{z) given by 
Eq. ([l2|). Since the expansion of Gu{z) is regular for z ^ 0, all the moments of 
Tu are finite. From the leading singularity of Gu{z) around z — ±i, one obtains 
the asymptotic behavior of !Fu{x) for large arguments as: 



Tu{x) 



-\x\ 



(13) 



In the special case where i/ = — 1, corresponding to a 'volatility random walk', 
the sum in ( |l2[ ) can be explicitly performed, and leads to: 

'sinh(7rz) ~ 



G_i(z)=log 



(14) 




Figure 3: Graph of the stable distribution P{x) = J-"_i(x). 
behavior !F_i{x) ~ exp(— |a;|) is shown as the dotted curve. 



The asymptotic 



This distribution is shown in Fig. 3, together with the predicted asymptotic 
behavior (dotted hne). The kurtosis of this distribution is equal to 6/5. Inter- 
estingly, J- 1, has a shape similar to hyperbolic distributions J23| with exponential 
tails which have been proposed in a financial context (see |18| ) . The appearance 
of such laws might thus be related to the existence of long-ranged correlations 
in the volatility. 

In summary, the purpose of this paper was to show, on an exactly soluble 
'stochastic volatility' model, that an apparent multiscaling behavior can appear 
as a result of very long transient effects, induced by the long range nature of 
the volatility correlations. This model is inspired by real price time series, and 
leads to an effective exponent spectrum (^q in close correspondence that reported 
in recent papers on the subject. We therefore suspect that indications of mut- 
lifractal behavior found in financial data might be misleading, as they could 
be caused by crossover effects that do not correspond to the true asymptotic 
behavior. To check more carefully for crossover effects, it might be helpful to 
analyze not only the moments but also the cumulants in empirical studies. We 
also have, en passant, found a new family of stable laws for sums of correlated 
random variables in the case where the volatility correlation is growing with 
time. It would be very interesting to characterize the attraction basin of these 
new stable laws. 
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